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Abstract:  The  theoretical  and  experimental  demonstration  of  a 
multispectral  Sagnac  interferometer  (MSI)  is  presented.  The  MSI  was 
created  by  including  two  multiple-order  blazed  diffraction  gratings  in  both 
arms  of  a  standard  polarization  Sagnac  interferometer  (PSI).  By  introducing 
these  high-order  diffractive  structures,  unique  spectral  passbands  can  be 
amplitude  modulated  onto  coincident  carrier  frequencies.  Extraction  of  the 
modulated  multispectral  images,  corresponding  to  each  passband,  is 
accomplished  within  the  Fourier  domain.  This  yields  a  unique  multispectral 
sensor  capable  of  imaging  all  the  passbands  in  a  single  snapshot.  First,  the 
theoretical  operating  principles  of  a  polarization  Sagnac  interferometer  are 
discussed  to  provide  a  context  for  the  MSI.  This  is  followed  by  the 
theoretical  and  experimental  development  of  the  MSI,  which  is  an  extension 
of  a  dispersion-compensated  polarization  Sagnac  interferometer  (DCPSI) 
based  polarimeter.  Indoor  and  outdoor  testing  and  validation  of  the  MSI  is 
performed  by  observing  vegetation,  demonstrating  the  ability  of  our 
experimental  setup  to  detect  four  distinctive  spectral  passbands. 
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1.  Introduction 

Multispectral  imaging  is  a  valuable  technique  in  remote  sensing  and  biomedical  imaging  [1, 
2].  Several  sensing  modalities  exist  for  measuring  multispectral  data.  One  method  time- 
sequentially  images  a  scene  through  multiple  filters,  placed  in  front  of  a  standard  lens  and 
focal  plane  array  (FPA)  [3].  This  produces  a  relatively  simple  and  inexpensive  instrument. 
However,  since  each  spectral  passband  is  measured  at  a  different  time,  error  is  introduced 
when  viewing  temporally  dynamic  scenes.  Consequently,  the  multispectral  passbands  must  be 
measured  in  parallel.  For  some  implementations,  these  parallel  measurements  are  obtained 
with  multiple  cameras,  lenses  and  filters  [4].  Another  method,  which  is  likely  the  most 
common  use  and  implementation  of  multispectral  imaging,  can  be  considered  a  division  of 
FPA  (DoF)  approach.  In  DoF,  different  filters  are  placed  onto  the  pixels  of  a  CCD  camera, 
similar  to  a  standard  consumer  color  camera  [5].  For  application-specific  selection  of  the 
passbands,  “super-pixels”  are  created  that  typically  contain  2x2  pixels,  where  each  pixel 
contains  a  unique  spectral  passband.  This  produces  a  compact  4  band  multispectral  camera 
[6]. 

The  method  for  multispectral  imaging,  demonstrated  here,  is  derived  from  a  snapshot 
imaging  polarimeter  initially  developed  by  K.  Oka  [7].  Oka’s  interferometric  technique 
demonstrated  that  the  2D  Stokes  parameters  can  be  amplitude  modulated  onto  spatial  carrier 
frequencies.  Reconstruction  of  the  2D  state  of  polarization  is  accomplished  by  Fourier 
filtering  each  carrier  frequency  component.  This  enables  snapshot  measurement  of  the 
complete  2D  Stokes  vector.  However,  one  drawback  of  this  technique  is  that  long  coherence- 
length  (narrow-band)  illumination  is  required.  For  applications  requiring  passive  illumination, 
limitations  are  encountered  due  to  the  low  signal  to  noise  ratio  of  the  measured  data  [8].  In  an 
attempt  to  remedy  the  narrow-bandwidth  limitation,  the  dispersion-compensated  polarization 
Sagnac  interferometer  (DCPSI)  was  developed  [9] .  Through  its  mathematical  development,  it 
was  discovered  that  a  further  extension  of  the  carrier  frequency  technique  can  be  realized. 
Whereas  the  original  DCPSI  amplitude  modulates  Stokes  parameters  onto  carrier  frequencies, 
it  is  possible,  with  a  minor  modification,  to  modulate  unique  spectral  passbands  onto  carrier 
frequencies.  This  device  creates  a  new  sensor  implementation  for  snapshot  multispectral 
imaging. 

In  this  paper,  we  outline  the  theoretical  and  experimental  development  of  a  multispectral 
Sagnac  interferometer  (MSI).  In  section  2,  the  theoretical  background  of  a  standard 
polarization  Sagnac  interferometer  (PSI),  when  used  for  polarimetric  imaging,  is  outlined  to 
provide  a  context  for  the  work.  In  section  3,  details  regarding  the  theory  of  the  DCPSI  are 
overviewed.  In  section  4,  the  multispectral  variant  of  the  DCPSI  is  introduced,  and  the  key 
differences  between  the  MSI  and  DCPSI  are  explained.  Lastly,  in  section  5,  the  experimental 
validation  and  calibration  of  the  MSI  are  detailed,  followed  by  a  brief  conclusion  of  the  work 
in  section  6. 


2.  Polarization  Sagnac  interferometer  for  polarimetric  imaging 

The  multispectral  Sagnac  interferometer  (MSI)  is  an  extension  of  a  standard  polarization 
Sagnac  interferometer  (PSI)  when  it  is  used  for  Stokes  imaging  polarimetry  [10].  One 
configuration  for  a  PSI  is  depicted  in  Fig.  1. 


Fig.  1.  One  configuration  for  a  polarization  Sagnac  interferometer.  The  distance  between  the 
WGBS  and  mirrors,  Mi  and  M2,  is  denoted  d\  and  d2,  respectively.  A  shear  (SPSi)  is  produced 
when  d\  ±  d2.  The  case  d\  >  d2  is  illustrated,  with  d\=d2  +  a. 

The  PSI  consists  of  two  mirrors,  Mi  and  M2,  with  a  wire  grid  polarization  beamsplitter 
(WGBS),  a  focal  plane  array  (FPA),  and  an  objective  lens  with  focal  length /obj.  The  WGBS  is 
oriented  with  its  transmission  axis  at  0°  with  respect  to  the  v-axis.  Two  beams,  with  a  shear 
Spsi,  are  created  when  d\  ±  d2.  If  the  distance,  d\,  is  offset  by  an  amount  a ,  such  that  d\  =  d2  + 
a ,  the  shear  is 


Sps,=d2a.  (1) 

Therefore,  the  shear  is  directly  proportional  to  the  difference  in  the  mirror  to  WGBS 
separation.  Furthermore,  the  shear  is  achromatic;  it  is  independent  of  the  wavelength  of  the 
incident  illumination. 

The  functional  form  of  the  intensity  pattern  incident  on  the  FPA  can  be  calculated  by 
considering  the  optical  path  difference  (OPD).  The  OPD  between  the  two  sheared  beams 
exiting  the  system  is  depicted  in  Fig.  2,  and  is  given  by 

OPD  =  SPSI  sin(#)  ~  SPS10  ,  (2) 


where  0  is  the  angle  of  the  incident  ray  [11].  When  the  two  beams  are  combined  by  the  lens, 
they  produce  interference  fringes  on  the  FPA.  The  interference  between  the  two  sheared  rays 
can  be  expressed  as 


I  PSI  {xi^yi)  ( 


Ex{xi,y,,t)e 
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where  <  >  represents  the  time  average,  xt  and  yt  are  image-plane  coordinates,  and  cph  (p2,  are 
the  cumulative  phases  of  each  ray.  Expansion  of  this  expression  yields 
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where  E*,  Fy  are  now  implicitly  dependent  on  xt  and  The  phase  factors  are 
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Fig.  2.  OPD  generated  by  a  shearing  distance,  SPSi,  for  a  PSI  and  simplified  SPP. 

Ipsi  can  be  re-expressed  in  terms  of  the  Stokes  parameters.  The  Stokes  parameters  are  defined 
from  the  components  of  the  electric  field  as  [12] 
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Re-expressing  Eq.  4  by  use  of  the  Stokes  parameter  definitions  and  (p h  (p2,  yields 


Ipsi  ixi^i )  ~~  2 
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(7) 


Consequently,  the  PSI  modulates  the  Stokes  parameters  S2  and  S3  onto  a  carrier  frequency, 
while  Sq  remains  as  an  un-modulated  component.  Fourier  filtering  can  then  be  utilized  to 
demodulate  and  reconstruct  the  spatially  dependent  Stokes  parameters.  The  spatial  carrier 
frequency  modulating  the  Stokes  parameters  is 


U  psi 


V2 a 

4 fobj  ' 


(8) 


Notable  in  the  carrier  frequency  is  the  inverse  relationship  to  the  wavelength.  Therefore, 
spectrally  narrow-band  illumination  must  be  used  in  order  to  maintain  high  carrier  frequency 
visibility  across  the  FPA.  This  narrow-bandwidth  requirement  can  be  disadvantageous  with 
regard  to  the  system’s  radiometric  throughput,  thereby  decreasing  the  signal-to-noise  ratio  of 
the  reconstructed  Stokes  parameters. 

To  maintain  high  fringe  visibility  over  a  broad  bandwidth  (white  light),  a  dispersion 
compensated  polarization  Sagnac  interferometer  (DCPSI)  was  developed.  The  DCPSI 
generates  a  shear  that  is  linearly  proportional  to  the  wavelength,  such  that  Sshear  Ay,  where 
y  is  an  optical  thickness,  analogous  to  a. 


3.  Dispersion  compensated  polarization  Sagnac  interferometer  (DCPSI) 

Compensation  of  the  dispersion  in  the  PSI’s  carrier  frequency  can  be  realized  with  the 
introduction  of  two  blazed  diffraction  gratings  [9].  Diffraction  gratings  are  known  for  their 
ability  to  generate  white-light  interference  fringes  [13].  The  optical  layout  now  takes  the  form 
of  the  dispersion-compensated  PSI  (DCPSI)  depicted  in  Fig.  3,  and  is  the  PSI  observed 
previously  in  Fig.  1  with  =  d2  and  the  inclusion  of  two  identical  gratings,  Gi  and  G2.  A 
ray’s  diffraction  angle,  after  transmission  through  Gi  or  G2,  is  calculated  for  normal  incidence 
by 


sin(<9)  = 


mA 
~d  ’ 


(9) 


where  6  is  the  diffraction  angle  of  the  ray  as  measured  from  the  grating’s  normal,  m  is  the 
order  of  diffraction,  and  d  is  the  period  of  the  grating.  Since  d  is  typically  large  (>  30  pm), 
small  angle  approximations  can  be  used  to  simplify  Eq.  9,  yielding  0  ~  mAfd  . 


Fig.  3.  DCPSI  with  blazed  diffraction  gratings,  Gi  and  G2,  positioned  at  each  output  of  the 
WGBS.  To  remove  the  achromatic  shear,  the  distance  between  the  WGBS  and  mirrors,  Mi  and 
M2,  depicted  previously  in  Fig.  1,  is  d\  =  d2.  Inclusion  of  the  gratings  generates  a  shear  that  is 
directly  proportional  to  the  wavelength. 

In  the  DCPSI,  the  spectrally  broadband  beam  transmitted  by  the  WGBS  is  diffracted  by  Gi 
into  the  +1  order.  When  these  dispersed  rays  are  incident  on  G2,  the  diffraction  angle,  induced 
previously  by  Gb  is  removed.  The  rays  emerge  parallel  to  the  optical  axis,  but  are  offset  by  a 
distance  proportional  to  -Ax0,  where  is  some  constant  related  to  the  DCPSI’ s  parameters. 
Conversely,  the  beam  reflected  by  the  WGBS  is  initially  diffracted  by  G2.  The  dispersed  rays 
are  then  diffracted  to  be  parallel  to  the  optical  axis  by  Gb  and  exit  the  system  offset  by  a 
distance  proportional  to  +Ax0.  Assuming  small  angles  and  that  Gi  and  G2  are  identical,  the 
shear,  S dcpsii  is 


Sdcpsi  =2—(a  +  b  +  c)  (10) 

d 

where  a ,  b ,  and  c  represent  the  distances  between  Gi  and  Mb  Mi  and  M2,  and  M2  and  G2, 
respectively.  Hence,  the  DCPSI  can,  to  first  order,  generate  a  shear  that  is  directly 
proportional  to  the  wavelength.  This  makes  the  phase  factors 
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Using  these  phase  factors,  the  intensity  on  the  FPA  is 
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where  the  total  intensity  pattern  is  a  summation  from  the  minimum  to  maximum  order  of  a 
blazed  grating,  such  that  the  maximum  order  is  (J/2min)sin(7i/2),  where  2min  is  the  minimum 
wavelength  passed  by  the  optical  system.  Sq  (m),  S2  (m),  and  S2  (m)  are  the  Stokes  parameters 
weighted  by  the  diffraction  efficiency  (DE)  of  both  gratings  before  integration  over 
wavelength, 


S0(m)=  J  DE2(A,m)S0(A)dA 

-U 

S2(m)  =  J  DE2(A,m)S2(A)dA 

^min 

imax 
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(15) 


where  2min  and  2max  denote  the  minimum  and  maximum  wavelengths  passed  by  the  optical 
system.  The  carrier  frequency,  UDCpsh  is 

2m 

UDCPSI=-^(a+b  +  c).  (16) 

dLj 

Therefore,  insertion  of  the  blazed  diffraction  gratings  enables  the  dispersion,  observed 
previously  in  Eq.  8,  to  be  removed  from  the  carrier  frequency.  Application  of  the  DCPSI 
towards  imaging  polarimetry  is  described  in  Ref.  [9]  in  which  a  single-order  blazed  diffraction 
grating  is  utilized.  However,  the  carrier  frequency’s  dependence  upon  the  diffraction  order,  m, 
can  be  exploited  for  multispectral  imaging.  This  is  accomplished  by  substituting  the  single¬ 
order  gratings  of  the  DCPSI  with  multiple-order  gratings. 

4.  Multispectral  variant  of  the  DCPSI 

The  MSI  is  nearly  identical  to  the  DCPSI  depicted  in  Fig.  3,  except  that  the  blazed  gratings, 
Gi  and  G2,  now  have  multiple-orders  over  the  sensor’s  operating  bandwidth.  The  difference 
between  the  gratings  is  depicted  in  Fig.  4,  where  the  multiple-order  structure  in  Fig.  4  (b)  has 
a  taller  grove  profile  than  its  single-order  counterpart,  per  Fig.  4  (a).  The  theoretical 
diffraction  efficiency  (DE)  for  an  ideal  blazed  grating  is  calculated  by, 

DE(A,m)  =  sine2  ^  j  (17) 


with 


OPD  =  h(nl-n2)  (18) 

where  h  is  the  grove  height  and  nx,  n2  are  the  indices  of  refraction  for  the  incident  and  blaze 
medium,  respectively  [14]. 


(a)  (b) 


Fig.  4.  (a)  Single  order  blazed  grating,  where  d  is  the  period,  n\  and  n2  are  the  indices  of 
refraction  for  the  incident  and  blaze  medium,  respectively,  and  hi  is  the  depth  ( OPD  ~  1  wave). 

(b)  Multiple-order  blaze  grating,  where  h2  is  typically  3-10  times  larger  than  h\. 

Using  Eq.  17,  the  theoretical  diffraction  efficiency  for  both  a  single  and  multiple  order  blazed 
grating  was  calculated.  Fig.  5  (a)  demonstrates  the  efficiency  assuming  a  single-order  grating 
with  hx  =  1.28  pm,  nx  -  1.0  and  n2  =  1.5.  This  yields  a  first  order  blaze  wavelength  of  Xx  = 
640  nm,  making  it  the  dominant  order  for  wavelengths  spanning  425  -  1000  nm.  Similarly, 
Fig.  5  (b)  depicts  the  efficiency  for  a  multiple-order  grating  with  h2  =  4.07  pm,  nx  =  1.0  and  n2 
=  1.5.  This  produces  a  first  order  blaze  wavelength  of  Xx  -  2035  nm,  meaning  that  the  higher 
orders  (m  =  2  through  5)  are  distributed  over  400-1000  nm.  The  blaze  wavelength  for  the 
higher  orders  is  Xx /m,  such  that  22  =  1017  nm,  23  =  678  nm,  24  =  508  nm,  and  25  =  407  nm. 


Wauslength  (nm)  Wavelength  (nm) 
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Fig.  5.  (a)  Diffraction  efficiency  (percent)  for  a  single  order  blazed  grating  in  air  with  h\  =  1.28 
pm  and  n  =  1.5.  (b)  Multiple  order  blazed  grating  in  air  with  h2  =  4.07  pm  and  n  =  1.5. 

As  is  demonstrated  in  Fig.  5  (b),  use  of  a  multiple-order  grating  will  generate  several 
diffraction  orders  within  the  system,  where  each  order  has  a  unique  blaze  wavelength  and 
spectral  passband.  Due  to  the  carrier  frequency’s  dependence  on  the  order  m,  each  spectral 
passband  is  modulated  onto  a  different  spatial  carrier  frequency.  Isolating  each  passband  is 
accomplished  by  use  of  Fourier  filtering,  detailed  in  section  5.1. 

Besides  including  the  multiple-order  gratings,  an  additional  modification  is  needed  to  the 
system  (Fig.  3)  for  measurements  of  the  incident  spectral  components;  namely,  a  linear 
polarizer  (FP)  must  be  inserted  in  front  of  the  WGBS.  If  the  FP’s  transmission  axis  is  oriented 
at  45°  with  respect  to  the  v-axis,  then  the  Stokes  vector  incident  on  the  WGBS  is 
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The  intensity  pattern  for  the  MSI  becomes 
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It  should  be  noted  that  the  dominant  orders  experimentally  observed  in  the  system 
correspond  to  the  ceiling  (Ce)  of  2i/2min,  where  X\  is  the  first  order  blaze  wavelength  of  the 
diffraction  grating.  This  changes  the  maximum  limit  of  the  summation  from  d/2m in  to 
Ce[2i/2min]  in  Eq.  20.  Furthermore,  if  S2  is  negligible  in  Eq.  21,  then  the  modulated 
information  corresponds  to  the  incident  50  component  (i.e.  the  total  incident  irradiance). 
Consequently,  S2  in  the  scene  can  cause  error  in  the  resulting  multispectral  data.  To  make  the 
contribution  of  S2  negligible,  scenes  consisted  of  diffuse  surfaces  oriented  normal  to  the 
optical  system.  Additionally,  the  illumination  source  was  always  oriented  to  preferentially 
reflect  light  in  the  S 1  state.  Further  analysis  of  the  error  due  to  this  effect  is  beyond  the  scope 
of  the  current  study.  However,  a  design  that  is  effectively  polarization  insensitive  can  be  made 
by  replacing  the  WGBS  with  a  standard  50/50  beamsplitter. 

5.  Experimental  verification  of  the  MSI 

To  verify  the  operating  principles  of  the  MSI,  the  experimental  setup  depicted  in  Fig.  6  was 
implemented.  Spatially  and  temporally  incoherent  light  is  configured  by  aiming  the  output  of 
a  fiber-light,  sourced  by  a  tungsten-halogen  lamp,  onto  a  diffuse  object.  A  dichroic  polymer 
linear  polarizer  at  45°  is  used  to  linearly  polarize  the  light,  such  that  the  So  Stokes  parameter 
can  be  modulated  onto  the  carrier  frequencies  (Eq.  19).  The  WGBS  has  a  clear  aperture  of  21 
mm  and  consists  of  anti-reflection  coated  aluminum  wires.  Mirrors  and  M2  are  25  mm 
diameter  l/10th  wave  optical  flats.  A  zoom  lens,  with  focal  length  fc  =  70-300  mm  (F/ll),  is 
used  to  image  the  object  plane  to  infinity,  while  an  objective  lens,  with  focal  length  fobj  =  250 
mm,  is  used  for  re-imaging.  The  diffraction  gratings,  Gi  and  G2,  are  blazed  for  a  first  order  ( m 
=  1)  wavelength  of  =  2035  nm  on  a  BK7  substrate  with  a  period  d  =  69.9  jum.  A  depiction 
of  the  ideal  theoretical  diffraction  efficiencies  for  these  grating  parameters  was  provided 
previously  in  Fig.  5  (b).  Fastly,  an  IR  blocking  filter  is  used  to  reject  light  with  wavelengths 
longer  than  750  nm. 

Another  aspect  of  the  setup  involves  the  placement  of  an  iris,  set  to  a  4  mm  diameter, 
between  the  zoom  lens  and  WGBS.  During  initial  experimentation,  it  was  observed  that  the 
fringe  visibility  decreased  as  the  aperture  diameter  increased.  Consequently,  the  optical 
surface  of  the  WGBS  was  quantified  in  a  T wyman-Green  interferometer  to  have  3  waves  of 
astigmatism.  Therefore,  the  iris  restricts  the  illuminated  area  of  the  WGBS  to  reduce  the  effect 
of  this  aberration. 


Fig.  6.  Experimental  setup  of  the  MSI  with  a,  b,  and  c  set  to  13.3  mm,  51.4  mm,  and  19.1  mm, 
respectively.  An  IR  blocking  filter  keeps  wavelengths  >  750  nm  from  entering  the  system. 

5. 1  MSI  calibration 

In  order  to  calibrate  and  reconstruct  data  from  the  instrument,  data  processing  is  accomplished 
in  the  Fourier  domain  [7].  Taking  the  Fourier  transform  of  the  intensity  pattern  (Eq.  12)  yields 
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where  3  is  the  Dirac-delta  function  and  £,  rj  are  the  Fourier  transform  variables  along  xt  and  yh 
respectively.  Extraction  of  the  Stokes  parameters  is  accomplished  by  taking  an  inverse  Fourier 
transform  of  the  filtered  carrier  frequencies,  or  “channels.”  Filtration  and  inverse  Fourier 

transformation  of  the  S(£,r/)  (channel  C0)  and  S^  +  UMSI  (m),r/)  (channel  C\(m)) 
components  yields 

1  Ce[V4m„] 

*r‘[c0]=-  z  soH  (24) 

^  m= 0 

&~l  [ Cl (»»)]  =Is'0(m)exp(j2xUMSI  (m)xi ) .  (25) 

Extraction  of  the  spectral  passbands,  Sq  (m),  requires  demodulation  of  the  exponential  phase 
factor.  Since  knowledge  of  the  phase  is  irrelevant  for  reconstruction  (50  is  always  positive), 
demodulation  is  achieved  by  taking  the  absolute  value  of  the  filtered  channel, 

|lr-I[c1(m)]|=Z^(m) . 


(26) 


A  150x150  pixel  section,  from  a  raw  image  of  a  uniformly  illuminated  high  reflectance 
diffuser,  is  depicted  in  Fig.  7  (a).  The  Fourier  transformation  of  the  image  is  depicted  in  Fig.  7 
(b),  demonstrating  the  presence  of  a  center-burst  and  four  carrier  frequencies,  along  with  their 
conjugates. 

Notable  is  that  the  image  and  Fourier  domain  is  rotated  by  approximately  31°.  This 
rotation  was  induced  by  rotating  the  camera  to  reduce  reconstruction  artifacts  caused  by 
aliasing  between  adjacent  channels  [15].  Specifically,  it  reduces  aliasing  caused  by  the  large 
side-lobes  of  the  carrier  frequency  components  along  and  rj ',  originating  from  the  nature  of 
the  2D  sine  function.  Further  discussion  of  error  related  to  this  phenomenon  is  beyond  the 
scope  of  the  current  paper.  However,  it  should  be  mentioned  that  while  changing  the  shearing 
direction  {i.e.  rotating  the  gratings)  will  rotate  the  fringes  without  rotating  the  image,  doing  so 
was  not  feasible  in  our  experimental  setup  due  to  optomechanical  limitations. 


Fig.  7.  (a)  Raw  image  of  a  uniformly  illuminated  diffuser  demonstrating  the  fringe  pattern  in 
the  MSI.  (b)  Fourier  transformation  of  the  raw  image  data  in  (a).  Four  carrier  frequencies  are 
observed,  corresponding  to  m  =  2,  3,  4  and  5. 

5.2  Spectral  passbands 

The  operating  principles  discussed  previously  indicate  that  the  order-dependent  carrier 
frequencies  are  visible  only  when  the  sensor  is  illuminated  with  light  at  specific  wavelengths. 
To  validate  that  the  passbands  correspond  to  a  given  diffraction  order  and  to  measure  their 
relative  spectral  responsivity,  the  experimental  setup  depicted  in  Fig.  8  was  implemented.  A 
monochrometer,  sourced  by  a  tungsten  halogen  lamp,  is  used  to  generate  narrowband  light. 
The  exit  slit  of  the  monochrometer  was  then  imaged  by  the  MSI.  Using  a  focal  length  of  70 
mm  for  the  collimating  lens  ensured  that  the  edges  of  the  slit  were  not  visible,  such  that 
fringes  were  present  across  the  entire  FOV. 


Fig.  8.  Monochrometer  configuration  for  sending  light  into  the  MSI  for  verification  of  the 
fringe  visibility.  The  bandwidth  of  the  light  exiting  the  monochrometer  was  approximately  10.5 
nm  using  a  3  mm  exit  slit. 

To  calculate  the  relative  responsivity,  a  calibrated  spectrometer  was  used  to  characterize  the 
output  of  the  monochrometer.  For  radiometric  calibration  of  the  spectrometer,  a  NIST- 


traceable  tungsten  halogen  lamp  was  utilized.  This  enabled  the  output  of  the  monochrometer 
to  be  characterized  in  radiometric  units,  as  illustrated  in  Fig.  9  (Left).  The  response  of  each 
channel  was  then  measured  by  Eq.  26  as  the  monochrometer  was  scanned  from  400-800  nm  in 
10  nm  increments.  An  average  of  approximately  100x100  pixels  from  the  center  of  the  image 
was  then  calculated  and  normalized  by  the  measured  response  of  the  monochrometer.  The 
relative  response  for  each  channel  is  depicted  in  Fig.  9  (Right)  alongside  the  combined 
relative  response  of  the  FPA  and  IR  blocking  filter;  the  response  of  the  other  optical  elements 
in  the  system  (dichroic  linear  polarizers,  lenses,  beamsplitter,  etc.)  is  not  included  in  this 
curve.  Despite  this,  the  relative  spectral  response  of  each  channel  generally  follows  the  trends 
of  the  FPA  and  filter. 

Monochometer  Output  vs  Wavelength  Relative  Spectral  Response  vs.  Wavelength 
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Fig.  9.  (Left)  Measured  monochrometer  output,  proportional  to  W/m2.  (Right)  Measured 
relative  spectral  response  of  the  different  passbands  within  the  MSI.  Also  included  (solid  light- 
gray  line)  is  the  relative  response  of  the  FPA  multiplied  by  the  measured  transmission  of  the  IR 
blocking  filter. 

5.3  Multispectral  reconstruction  validation 

To  validate  the  accuracy  of  the  multispectral  reconstructions  from  the  MSI,  experiments  were 
performed  that  compared  MSI  results  to  that  of  the  calibrated  spectrometer  (U2S). 
Specifically,  reflectance  measurements  of  a  healthy  and  unhealthy  tree  leaf,  taken  from  an 
Acacia  crassifolia ,  were  measured  with  the  MSI  and  the  U2S.  These  samples  enabled  the 
characteristic  reflectance  features  of  chlorophyll,  which  is  present  in  the  healthy  leaf  but 
largely  absent  in  the  unhealthy  leaf,  to  be  used  for  the  comparison  [16,  17].  Reflectance 
measurements  from  the  U2S  were  matched  to  the  MSI’s  passbands  by 
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where  2  is  the  wavelength  in  nm,  TMSI  is  the  relative  response  of  the  passband  of  order  m  per 
Fig.  9,  and  Iieaf ,  I  diffuser  is  the  measured  irradiance  of  the  leaf  and  of  a  white  diffuse  reflectance 
panel,  respectively.  The  reflectance  measurement  from  the  MSI  is  calculated  by 
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where  /  and  n  are  integers  representing  pixels  in  x  and  y,  respectively.  Here,  the  inner  90%  of 
the  image  is  used  in  the  summation  to  average  over  a  large  area  (~25  mm2)  of  the  sample. 


The  two  experimental  configurations  for  the  reflectance  measurements  are  depicted  in 
Fig.  10,  where  the  source  and  viewing  geometry  between  both  the  U2S  and  MSI  is  maintained 
for  both  measurements.  Since  the  U2S  is  non-imaging,  defocus  was  implemented  in  the  MSI 
to  obtain  reflectance  measurements  that  average  over  the  observed  area  of  the  sample. 


Fig.  10.  (Left)  System  setup  for  U2S  reflectance  measurements.  The  incidence  angle  of  the 
source  (tungsten-halogen  lamp)  was  approximately  37°  with  respect  to  the  surface  normal  of 
the  sample.  (Right)  System  setup  for  MSI  reflectance  measurements.  The  sample  is  defocused 
to  average  the  reflectance  over  the  observed  area. 

The  measured  reflectance  from  the  U2S  and  the  MSI  yields  the  results  portrayed  in  Fig.  11. 


Relative  Reflectance  vs.  Order  (Band) 


Fig.  11.  Relative  reflectance  of  the  healthy  and  unhealthy  leaves,  normalized  to  the  irradiance 
in  order  2,  as  measured  with  the  MSI  and  U2S.  The  healthy  vegetation  experiences  more 
absorption  in  order  3  due  to  the  presence  of  chlorophyll. 

The  root  mean  square  (RMS)  error  between  the  MSI  and  U2S  data  is  calculated  with 

f=\lZ(*U2S  W-^MS/W)2-  (29) 

V  4  m= 2 


The  RMS  error  for  the  healthy  and  unhealthy  leaf  samples  yields  0.0115  and  0.0117, 
respectively.  Consequently,  the  error  in  the  spectral  reconstruction  is  on  the  order  of  1.2 
percent  between  the  U2S  and  the  MSI. 

5.4  Multispectral  spatial  reconstructions 

To  demonstrate  the  in  focus  imaging  capabilities  of  the  MSI,  both  the  healthy  and  unhealthy 
leaf  samples  were  imaged  simultaneously  with  the  sensor.  A  relative  reflectance  calculation 
was  performed  using 


RMSI{l,n,m)  = 


1 diffuser  ‘ 


(30) 


The  band-integrated  image,  taken  by  demodulating  the  first  component  in  Eq.  22,  is  depicted 
in  Fig.  12  (a),  while  an  image  from  each  band  (or  order)  is  depicted  in  the  bottom  row.  An 


additional  calculation,  demonstrating  the  applicability  of  the  results,  is  the  normalized 
difference  vegetation  index  (NDVI)  provided  in  Fig.  12  (b).  The  NDVI  is  calculated  by 


NDVI  (l,n) 


I(l,n,2)-I(l,n,3) 
l(l,n,2)  +  l{l,n,3)' 


(31) 


If  the  leaf  has  a  high  chlorophyll  content,  order  3  will  contain  less  reflected  energy  than  order 
2,  yielding  a  high  NDVI.  Conversely,  if  the  leaf  has  little  to  no  chlorophyll,  order  3  will  have 
nearly  equal  reflected  energy  compared  to  order  2,  yielding  a  low  NDVI.  This  is  observed  in 
the  NDVI  image,  where  the  upper  left  portion  of  the  scene  (quadrant  2)  contains  the  unhealthy 
leaf,  while  the  lower  right  region  (quadrant  4)  contains  the  healthy  leaf. 


m  =  0  NDVI 


Fig.  12.  Relative  reflectance  images  of  a  healthy  and  unhealthy  leaf,  (a)  Band-integrated  image, 
(b)  normalized  difference  vegetation  index  (NDVI).  (Bottom)  Image  from  each  order  (m  =  2 
through  m  =  5). 


Since  the  near  infrared  (NIR)  band  (order  2)  spans  wavelengths  from  685-755  nm  full  width 
half-max,  order  2  integrates  across  the  vegetative  red  edge.  Consequently,  the  NDVI 
calculation  performed  here  is  more  analogous  to  the  calculation  performed  in  Ref.  [17]. 
Gitelson  demonstrates  that  an  NDVI  calculation  performed  between  two  bands,  located  at 
525-555  nm  and  695-705  nm,  can  be  linearly  related  to  the  vegetation’s  chlorophyll 
concentration.  This  is  in  opposition  to  an  NDVI  calculation  performed  with  Landsat  7’s  band 
4  (780-900  nm)  with  ETM+  [18].  The  MSI  and  Landsat  7  NIR  band  is  overlaid  onto  the 
measured  relative  reflectance  for  Acacia  crassifolia  per  Fig.  13.  Ultimately,  the  system  can  be 
designed  to  incorporate  the  appropriate  bands  to  perform  an  NDVI  measurement  analogous  to 
Landsat  7,  or  other  passbands  given  the  requirements  of  a  specific  application. 

Application  of  the  MSI’s  multispectral  imaging  capabilities  is  further  illustrated  in 
outdoor  test  results.  A  photo  of  an  outdoor  scene,  taken  with  a  color  digital  camera,  is 
depicted  in  Fig.  14.  The  reconstructed  spectrally  band-integrated  image  from  the  MSI  is 
depicted  in  Fig.  15  (a),  while  an  image  from  each  band  is  depicted  in  the  bottom  row.  The 
NDVI  image  is  portrayed  in  Fig.  15  (b).  This  demonstrates  the  sensor’s  ability  to  operate  in  a 
remote  sensing  capacity. 


Relative  Reflectance  vs.  Wavelength  (Acacia  crassifolia) 
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Fig.  13.  Relative  reflectance  of  a  healthy  leaf  from  Acacia  crassifolia  measured  with  the  U2S 
(solid  black  line).  The  “red  edge”  begins  at  700  nm  and  peaks  at  approximately  765  nm.  The 
Landsat  7  (band  4)  spans  775-900  nm  (dark  dashed  line)  while  the  MSI  spans  685-755  nm 
(gray  dotted  line). 


Fig.  14.  Photo  taken  with  a  standard  color  digital  camera  of  an  outdoor  scene.  Healthy 
vegetation  is  present,  in  addition  to  brick,  concrete,  and  a  relatively  clear  sky. 


m  =  0  NDVI 


Fig.  15.  Relative  reflectance  images  of  an  outdoor  scene,  (a)  Band-integrated  image,  (b) 
normalized  difference  vegetation  index  (NDVI).  (Bottom)  Image  from  each  order  (m  =  2 
through  m  =  5). 


6.  Conclusion 


The  development,  calibration,  validation  and  experimental  results  of  a  different  kind  of 
multispectral  imager  have  been  presented.  The  multispectral  Sagnac  interferometer  (MSI)  was 
created  by  modifying  the  dispersion  compensated  polarization  Sagnac  interferometer  (DCPSI) 
originally  used  for  snapshot  polarimetric  imaging.  By  replacing  the  single-order  blazed 
diffraction  gratings  with  multiple-order  gratings,  each  diffraction  order’s  distinctive  spectral 
response  is  encoded  onto  unique  spatial  carrier  frequencies.  The  multispectral  reconstruction’s 
RMS  error,  as  calculated  between  the  MSI  and  a  calibrated  spectrometer,  was  measured  to  be 
approximately  1.2  percent  for  both  a  healthy  and  unhealthy  leaf  from  Acacia  crassifolia.  The 
imaging  potential  of  the  MSI  was  demonstrated  through  NDVI  measurements  of  the  same  leaf 
samples,  in  addition  to  outdoor  measurements  of  scenes  containing  vegetation.  Ultimately,  the 
sensor  compromises  spatial  resolution  to  capture  the  spectral  components.  However,  it  is 
envisioned  that,  with  further  optimization,  the  spatial  and  spectral  resolution  can  be 
maximized. 


